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A summary is given of the present status of the theory and experiment of the anomalous magnetic 
moment of the muon. A difference between predicted and measured values is an indication of physics 
beyond the Standard Model. A new experimental measurement has produced a value which differs 
from a recent Standard Model prediction by about 1.6 standard deviations. When first announced, 
the discrepancy was about 2.6 standard deviations, but theorists have recently found an error in the 
sign of the largest term in the standard model hadronic light-by-light contribution which reduces the 
difference. Additional data are being analyzed and elements of the theory are being scrutinized to 
provide, in the future, a sharper test of theory. 



1 Introduction 

The magnetic moment of a particle is given 
in terms of its spin by 

2mc 

and the anomaly is defined as a = ^{g — 2). 

Historically, the measurement of particle 
magnetic moments has been a valuable test 
of existing theories. For example, magnetic 
moment measurements on the hyperons have 
provided essential information on their sub- 
structure, and the electron anomaly has been 
the most stringent test of QED. 

The Dirac theory predicts that g = 2 
(a = 0) for point particles with spin ^ . While 
the hyperons have g factors very different 
from 2 because of their complex substructure, 
the leptons have g ~ 2 and anomalies which 
are nearly zero, consistent with the current 
evidence that they are point particles. The 
Standard Model predicts lepton anomalies on 
the order of one part in 800 due to their field 
interactions. In the cases of the electron and 
the muon, both the Standard Model predic- 
tions and the measurements are extremely 
precise, a relatively rare situation resulting 
in valuable tests of the theory. As we shall 



see, however, the muon typically has a much 
stronger sensitivity than the electron to any 
physics which has not been included in the 
Standard Model. 

By far the largest contribution to the 
lepton anomalies comes from the lowest or- 
der electromagnetic diagram, the Schwinger 
term (left diagram in Fig. ^), which gives 
a{QED; 1) = e.g. the same for muons 
and electrons. The next order electromag- 
netic diagrams, which involve virtual lepton 
(Fig. |l|) or hadron (Fig. |^) loops, are small 
compared to the Schwinger term, however, 
they are much larger for the muon than the 
electron as a result of the additional available 
rest mass energy. The difference in the con- 
tribution between the electron and muon in 
diagrams involving massive virtual particles 
typically scales as (^)^ ~ 40000, and it is 
this large factor which makes far more sen- 
sitive than to any unknown massive parti- 
cles. 

It is instructive to compare the values of 
the electrijn and muon anomalies. From Pen- 
ning trapa experiments, we have for the elec- 
tron al7 = (1159652.1884 ± .0043) x 10"^ 
(4 ppb) and positron a^^^ = (1159652. 1879± 
.0043) X 10^^ (4 ppb), agreeing within errors 
as required by CPT invariance. These val- 
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Figure 1. QED contributions to the anomaly. The 
first diagram is the lowest-order Schwinger term. 
The other diagrams are representative of higher-order 
QED contributions. 




Figure 2. The first-order hadronic diagram. 



ues are are in good agreement with the the- 
oretical prediction of QED to fourth order in 
f , ai? = (1159652.1535 ± .0240) x 10"^ (21 
ppbp. The theoretical error is dominated 
by the error in the value of a taken from 
Quantum Hall Effect experiments. If one as- 
sumes the correctness of QED, then the best 
determination of a comes from the electron 
anomaly measurements. The hadronic and 
electroweak contributions, 1.63(3) x 10~^^ 
and 0.030 x 10~^^, respectively, are small and 
the errors are negligible compared to the ex- 
perimental errors. 

The theoretical and experimental values 
for the are not known nearly as well as 
for Oe. In 1999, prior to the new published 
result, the world average of measured val- 
ues was af,'^ 
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Figure 3. Experimental measurements of a^. 



experiment (E821) usinff positive muon data 
taken in 1997 (13 ppm)i and 1998 (5 ppm)i. 
All of these measurements are in agreement 
within their errors, as shown in Fig. ^; the 
negative muon data were incorporated un- 
der the assumption of CPT invariance. Using 
the recent published compilation of the the- 
oretical irigrcdients to a^'* by Czarnecki and 
Marciano,! a^'' = (1165915.96 ± 0.67) x 10"^ 
(0.6 ppm), which used the hadronic evalu- 
ation by Davier and Hockerp we find that 
agreement between experiment and data was 
good: a^^P - a]^ = (4.5 ± 4.7) x 10'^. Re- 



Ironic light-by-light 
I which changes the 



cently, an error in the ha 
contribution was found,L 
theoretical prediction to a*^ — (1165917.68 ± 
0.67) X 10~^ (0.6 ppm). This improves the 
agreement between theory and experiment: 



-,th 



= (2.8 ±4.7) X 10" 



The new experimental result for 0^+ , 
based on a sample of 1 billion positrons col- 
lected in 1999 by E821,i is al='P{E82l) = 
(1165920.2±1.4±0.6) x 10"^ (1.3 ppm). The 



(1165920.5 ± 4.6) x 10 (4 new world value changes very little, ajf 



ppm). This included 9.4 and 10 ppm results 
for the positive and negative muon, respec- 
tively, in a series of famous experiments at 
CERN ending in the 1970's,i as weh as the 
results from the new Brookhaven muon (g-2) 



(1165920.3 ± 1.5) X lO'^ (1.3 ppm). Com- 
paring with the uncorrected theoretical num- 
ber gives a 2.6 a difference between measure- 



ment and theory: Aa^ — 



-,th 



= (4.3 ± 



1.6) X lO-y (3.7 ± 1.4 ppm). Including the 



light-by-light correction, Aa^^ = a^^^ — = 
(2.6 ± 1.6) X 10-9 (2.2 ± 1.4 ppm), a 1.6cr 
difference. 

Two additional data sets, from 2000 and 
2001 runs, are currently being analyzed. The 
2000 data set consists of about 4 billion 
events for the /i+ and 2001 consists of about 3 
billion events for the ii~ . The original stated 
goal of E821 was to reduce the anomaly mea- 
surement error to 0.35 ppm. Ideally, one 
would also equalize the errors on the /i"*" and 
fjL~ in order to optimally test CPT invariance 
and to study systematic issues in the experi- 
ment. We expect to come close to this goal, 
but this will require a future data run with 6 
billion events. 

2 Theory Status 

The Standard Model contributions to can 
be conveniently separated into QED, elec- 
troweak, and hadronic portions. Although 
the anomaly is dominated by the QED con- 
tribution, there are significant hadronic and 
EW contributions at the level of « 58.3 ppm 
and sa 1.3 ppm, respectively. We discuss 
each of these contributions below, with em- 
phasis on the hadronic contribution, whose 
error dominates the overall error in a*'*. 

The QED contribution, using a from 
al^P and calculated to fifth order in — (some 
of the highest order diagrams were esti- 
mated), contributesO an error of 21 ppb 
to af, = 116584705.7(2.9) x lO"". 

Kinoshitall3 reports that some of the fourth 
order terms are being re-calculated with 128 
bit precision, which may result in a small shift 
in the QED contribution; this is expected to 
have minimal impact on the comparison be- 
tween experimental and theoretical values of 

Op. 

The lowest order electroweak diagrams, 
involving the exchange of a W, Z or Higgs, 
are shown in Fig. ^. The electroweak con- 
tribution, including the 25% reduction from 
higher order terms, isi = 151(4) x 
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Figure 4. Lowest-order electroweak contributions. 



10-"(1.30 ± 0.03)ppm. The theoretical un- 
certainty is very small. The central value is 
right on the edge of the current experimental 
error, and it will be a significant contribution 
at the experimental error goal of 0.35 ppm. 

The lowest order hadronic diagram is 
shown in Fig. ^, where a hadron loop has 
been inserted into the Schwinger diagram. 
Since these contributions involve the strong 
interaction at low energies, they cannot be 
calculated from first principles. Their con- 
tribution can, however, be determined from 
measured e'^e" scattering cross sections over 
all energies through the use of the dispersion 
relation, Eq. ^ (also see Fig. ^(a)) 

a,ihad;l)^i^r T ^Kis)Ris) (1) 

K{s) is a slowly varying function and the 
e+e" data enter through the ratio of cross 

7->/ \ aie'^e^ — >hadrons) rr^i i 

sections, His) = , + , _/ . ihe low 

energy data are the most important as a re- 
sult of their large amplitude and the p- term 
in the integrand, where is the center-of- 
mass energy. 

High quality hadronic r decay data from 
LEP and Cornell can be used to augment the 
isovector part of the e^e^ data at energies 
below rrLrC^ (second diagram in Fig. ^) us- 
ing isospin invariance (to relate for example 
tt-tt" channels in r decays to tt'^t:^ channels 
in e+e" collisions) and the CVC hypothesis 
(to connect the W~ and photon intermediate 
states) .0 

There have been a number of evaluations 
of ap(had;l) over the past two decades. We 
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Figure 5. The electron scattering and tau decay di- 
agrams relevant to the determination of the lowest 
order contribution to a„. 



Table 1. As of 1995, contributions to ap(had;l) as a 
function of e+e~ energy, iikistrating the (soon-to-be 
obsolete) sources of error .E3. With the incorporation 
of new e+e~ and t data, the errors below 5 GeV 
decrease significantly from the values shown in the 
table. 





a^{had; 1) 


Error, ppm 


< 1.4 


87. % 


1.29 


1.4 -> 2.0 


4.6 % 


0.21 


2.0^ 3.1 


4.0 % 


0.30 


2.0 2.6 


2.9 % 


0.27 


2.6 ^ 3.1 


1.1 % 


0.12 


J/* (6 states) 


1.3 % 


0.08 


QCD 3.1 ^ oo 


3.0 % 


0.03 


Total 




1.37 



will look in some detail at those performed 



since 1995. The contributions of 



data 



to a^i{had; 1) in different energy ranges Jrom 
the evaluation by Brown and Worstell|13 are 
indicated in Table |[ (We note that the new 
data which are becoming available will render 
this table obsolete.) The largest contribution 
to the value and error are for ^/s < 1.4 GeV, 
a region dominated by the effects of the p 
resonance (see Figs. 6-8). The other energy 
ranges give a considerably smaller contribu- 
tion to the error, however for an 0.35 ppm 
measurement they cannot be neglected, par- 
ticularly in the range 1.4 2.6 GeV. 

Since 1995, there has been a substantial 
improvement in the data quality, but major 
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Figure 6. Preliminary pion form factor dEtta from 
Novosibirsk in the vicinity of the p resonance. t3 Inset: 
data the lu interference region. 



portions of the data are either preliminary or 
are still in the process of being incorporated 
into evaluations of a^t(had;l). 

Data on the pion form factor (which can 
be directly related to the e+e~ cross sec- 
tion) from the CMD2 and SND experiments 
at the VEPP-2M accelerator in Novosibirsk 
are nearing publication. They cover the im- 
portant energy range ^/s < 1.4 GeV. Prelim- 
inary data in the p resonance range (600 ta 
930 MeV) from CMD2 are shown in Fig. |.@ 
Their anticipated systematic error of 0.6% in 
this range would reduce the error contribu- 
tion from the e~^e~ data in this energy region 
by better than a factor of two. The VEPP- 
2000 project is an upgrade under construc- 
tion at Novosibirsk which will extend quality 
e~^e~ measurements up to 2 GeV, with an 
order of magnitude or more improvement in 
luminosity. 

Data relevant to the low energv region 
are also being taken by the KLOEo exper- 
iment at DA$NE. They operate at the </> 
mass, and then derive e^e~ cross sections 
at lower energies using the so-called radiative 
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Figure 7. Preliminary pion form facton-data from ra- 
diative return measurements (KLOE).Ej 



return method. In this approach, the initial 
state electron or positron radiates away some 
of its energy via a photon, providing access 
to scattering at energies below the (j) mass. 
Preliminary data are shown in Fig. |^. Cur- 
rently their systematic errors are on the or- 
der of a few percent, however, they expect to 
improve this in the near future to the point 
that they are competitive with the Novosi- 
birsk data. This measurement, along with 
possible plans at other e+e^ machines such 
as the B-factories and Cornell which have ac- 
cess to higher energies, will contribute signif- 
icantly to the next wave of precision e+e~ 
data. 

There are new, preliminary and pub- 
lished data from BESIItJ in the important 
energy range 2 — 5 GeV (Fig. ^). Note espe- 
cially that in the range 2 — 3 GeV, which has 
the largest contribution to a^(had;l), BESII 
data have much smaller errors and the cen- 
tral values are 15% lower compared to the old 
MARK I and Gamma2 data. 

High quality r decay data, which have 
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Figure 8. Preliminary and published e+e~ data, 2-5 
GeV, from BESII, along with older data from other 
experiment^lj • 



been produced over the last few years as de- 
scribed by P. Roudeau at this conference, 
have significantly reduced the contribution to 
the error of a^(had;l) for energies below the 
T mass.0 In addition, new analysis results 



should be avai' 
Some authorsllZI' 



ble soon from ALEPH.B 
J have questioned, however, 
whether there are sufficient controls over ra- 
diative corrections and corrections to the ap- 
proximations of CVC and isospin invariances. 
In addition, there appears to be an over- 
all normalization disagreement between the 
r data from ALEPH and Cornell.0 The r 
data represent, nevertheless, valuable addi- 
tions and checks on the e+e" data at low en- 
ergies, and various groups plan to study these 
issues in the future. 

A list of a^(had;l) evaluations is given 
in Table |^. Various combinations of ingredi- 
ents were used: e+e^ data from a wide range 
of experiments, r decay data from LEP and 
CESR, or theoretical input from perturbative 
QCD (pQCD) in the higher energy regions 
where it can be relied upon but where the 
data quality is poor. 

Eidelman and Jeger lehner0, EJ95, rehed 
primarily on e^e~ data, using pQCD only at 
the highest energies (> 20 GeV). Brown and 
WorstellEl, BW96, used essentially the same 
data set as EJ95, but took into account the 
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correlation of errors among data points com- 
ing from the same experiment. BW96 and 
EJ95 values and errors are in excellent agree- 
ment, suggesting that the correlation issue 
was not so important. Adel and YndurainEd, 
AY95, used pQCD in regions where the e+e" 
data were poor, and obtained a value some- 
what higher than, but still in agreement with, 
EJ95 and BW96. The theoretical input en- 
abled hem to quote a smaller £rror. In 1998, 
Adelman, Davier and H6cker,llil ADH98, did 
a full re-evaluation using updated values of 
the data sets used by EJ95 and BW96, ob- 
taining the same error and a central value 
which was lower but still consistent within 
errors. ADH98 then incorporated the high 
quality r data from ALEPH, producing a 
dramatic improvement in the error of the con- 
tribution below the r mass, and reducing the 
overall error in by 40%; inclusion of r data 
increased the central value somewhat. Subse- 
quently, Davier and H6ckerc3, DH98a, refined 
the e~^e~ and t evaluation of ADH98 by ap- 
plying pQCD above 1.8 GeV in regions where 
the data were poor, resulting in another 20% 
decrease in the error and a decrease in the 
central value. Most of the changes can be 
attributed to the effect of using pQCD from 
1.8 to 3 GeV, where the old data are rather 
poor, as seen in Fig. ||. The DH98a evalu- 
ation was done before the new BESII data 
(Fig. H) were available in this energy range; 
their pQCD calculations are in very good 
agreement with the BESII data, and about 
15% below the old data, pointing to the reli- 
ability of pQCD at these energies. In a sub- 
sequent work (DII98b) the same authors ap- 
plied QCD sum rules at low energies, result- 
ing in slight further reductions in error and 
central value. The value from DII98b was 
used in the Czarnecki and MarcianpB theo- 
retical compilation and also by E82llj to com- 
pare with their new experimental number. 

Some preliminary and published theoret- 
ical evaluations of a^(had;l) have appeared 
since the E821 publication of Brown, et al.i 



Table 2. A list of a number of recent evaluations of 
afj_{had\ 1). The last entry is the E821 experimental 
'measurement' of a^(had;l), obtained by subtracting 
the QED, EW and higher order hadronic contribu- 
tions from the experimental number. This is-^ up- 
dated version of a similar table found in RcfHfl. 



Ref. 


a^i{had] 1) 
(xlO") 


Comment 


EJ95N 


7024(153) 




BW960 


7026(160) 




AY95y 


7113(103) 


e+e-,QCD 


ADH980 


6950(150) 


e+e^ 


ADH980 


7011(94) 


e^e" , T 


DH98ay 


6951(75) 


pQCD 


DH98bl] 


6924(62) 


e+e^, r 
pQCD, sum rules 


NO 10 


7021(76) 


e^e~ , T 


TYO10 


6966(73) 


e'^e^ ,T 
space-like F^r 


EOliI 


6932(65) 




E82li 


7350(153) 


Expt -[ QED 
+EW-H(Had>l)] 
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NarisonEj, NOl, has used essentially the same 
T and e+e~ data sets at low energies and for 
the resonances as ADH98, with QCD applied 
to the continuum at the higher energies (> 
1.7 GeV), arriving at nearly the same value 
as ADH98, but with a slight y smaller er- 
ror. Troconiz and Yndurain§f, TYOl, have 
applied the maximum available data (includ- 
ing some preliminary data) and theory, fol- 
lowing the earlier approach of AY95, and also 
incorporated pion form factor data from pion 
scattering at low energies, to arrive at a value 
which compares closely to DH98b, with a 
slightly larger centraljialue and error. 

Finally, EidelmaiJl3, EOl, has produced a 
preliminary number based entirely on e^e~ 
data (except for pQCD at the very highest 
energies), including the new preliminary re- 
sults from Novosibirsk and BESII. He obtains 
a value and error which arc almost the same 
as DH98b. 

What can we conclude from this series 
of a^{had; 1) determinations? The first thing 
we note is that within their stated errors, all 
of the evaluations are in agreement. In par- 
ticular, the most recent evaluations are in ex- 
cellent agreement even with their smaller er- 
rors, although we note that the DH98b eval- 
uation has the smallest value leading to the 
largest discrepancy between theory and ex- 
periment. Secondly, an analysis which in- 
corporates all of the new r and e~^e~ data 
would be helpful. Eidelman, Davier and 
Hocker are presently collaborating on such an 
evaluation. When the final analyses, which 
include the latest excellent data are com- 
pleted, we can expect a more reliable value 
for a^{had; 1) with a smaller error. 

The higher order hadronic contributions, 
a^{had;> 1), can be separated into two 
parts. One part, involving higher order 
diagrams such as those in Fig. ^, have 
a relatively small contribution to and 
the error is negligible :EZl a^i{had]> 1) = 
— 101(6) X 10^^^. The other part, involv- 
ing the hadronic light-on-light diagram in 



Fig. |l^ (LOL), presents special problems be- 
cause, unlike other hadronic terms, it can- 
not be estimated based on experimental data^ 
The value used in Czarnecki and MarcianoEl 
is an average of the values aiid the errors of 
two separate determinationsc3' (which are 
in agreement within errors), a^(had;LOL) = 
—85(25) X 10~^^. Both calculations use mod- 
els motivated by chiral perturbation theory 
to calculate the contributions at low ener- 
gies. The largest contribution comes from 
the TT*^ pole term, with lesser contributions 
from the 7] and r/ poles. Other contributions 
can be w 20 — 30% in size relative to the 
pole terms, but when added together they 
largely cancel. After the publication of the 
new E821 experimental resultu and after this 
talk was given at LPOl, Knecht and Nyffeler 
E3 calculated the pion pole contribution us- 
ing Large-iVc and short-distance properties 
of QCD. They obtained virtually the_same 
magnitude as in references and but 
with the opposite sign: a|i/iad; LOL, tt") = 
+5.8(1.0) X IQ-i". The signs of the r] and ry' 
poles also change, however, Knecht and Nyf- 
feler only estimated their magnitude using a 
VMD model; their total pseudo-scalar con- 
tribution is a^{had; LOL, PS) = +8.3(1.2) x 
lQ-w_ The sign error hasLaince b£fin acknowl- 
edged iry the authors of Ej and Ej (see refer- 
ences u and E3). Changing the signs of the 
pole and axial vector terms, but keeping the 
other terms as calculated in l2I and Ej the 
same, increases the theoretical value of by 
17.2 X 10^^°, very close to just reversing the 
sign on the entire LOL contribution. Knecht 
and Nyffeler plan in the future to calculate 
all of the other terms in the hadronic LOL 
contribution. Several other groups are also 
considering new ways to tackle this difficult 
calculation. Ultimately, the hadronic LOL 
term may prove to be the limiting factor in 
the theoretical error of a^. 

The theoretical error on ap(/iac?;l) has 
gone down dramatically over the years, and 
with new, more accurate calculations of the 
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Figure 9. Some diagrams of the higher-order 
hadronic contributions. 




Figure 10. Diagram for the hadronic light-on-hght 
contribution. 



light-by-light contribution, the error and reli- 
ability of a^'' should continue to improve sig- 
nificantly. 

New physics will be reflected by a non- 

a' 



zero value of Aa"*^™ 



'^^P - a*'\ Some 



examples are muon substructure, anoma- 
lous gauge couplings, leptoquarks, or super- 
symmetry. In a minimal supersymmetric 
model with degenerate sparticle masses (see 
Fig. |ll|) , the contribution to Aa^^™ would be 
substantial in the case where tan /3 is large: 
^aSUSY ^ X l0-ii(Mf£li)2tan/3, 

where m is the sparticle mass. For 4 < 
tan/3 < 40, m w 150 - 500 GeV. Note 
that the supersymmetric diagrams (Fig ^) 
are analogous to the electroweak diagrams 




Figure 11. Lowest order diagrams for SUSY contri- 
butions. 



(Fig ^) . If supersymmetry is to explain all of 
Aa^*^™, then its contribution is large: almost 
two times bigger than the electroweak con- 
tribution. Of course, the presently observed 
Aa^ may also be due to a statistical varia- 
tion in the experimental number or to errors 
in the experiment or theory. When the addi- 
tional data are analyzed, and with the contin- 
ued extensive studies of the theory, one can 
anticipate that major progress will be made 
in understanding any non-zero Aa^. 

3 Experiment 

The ongoing muon (g-2) experiment at 
Brookhaven National Laboratory, E821, had 
its beginnings in the early 1980's. Its origi- 
nal goal was to measure to 0.35 ppm, or 
about 20 times better than the CERnH ex- 
periment. It received Laboratory approval in 
1987, and major construction on the storage 
ring magnet began in the early 1990's. The 
first data were taken in 1997, with one major 
run in each year 1998-2001. The experimen- 
tal technique follows the general one used in 
the CERN experiment with a number of im- 
portant improvements and innovations. 

The Brookhaven AGS delivers up to 7 x 
10^^ protons per bunch, with energies of 24 
GeV, onto a water-cooled, rotating nickel tar- 
get. There are 6-12 bunches per « 2.5 second 
AGS cycle, each about 50 ns wide FWHM 
and spaced 33 milliseconds apart. Secondary 
pions emitted from the target with momenta 
of 3.1 GeV/c are sent down a 72 m straight 
section of alternating magnetic quadrupoles, 
where highly polarized muons from forward 
pion decays are collected. The beam is then 
momentum-selected for either pions or the 
slightly lower- momentum muons, and then is 
injected through a field-free infiector@ region 
into a circular storage ring possessing a very 
homogeneous magnetic field. 

For the case of pion injection, with the 
pion momentum slightly higher than that of 
the central storage ring momentum, a small 



fraction 25 ppm) of the muons from pion 
decays will have the correct monienta and 
directions to be stored. The efficiency of 
this process is low, and the very high inten- 
sity of pions and secondary particles associ- 
ated with pion interactions with surround- 
ing materials creates severe background prob- 
lems for the detector system near the injec- 
tion time ('flash'). An essential improvement 
over prior experiments was the incorporation 
of direct muon injection. With muon in- 
jection, the number of stored muons is in- 
creased by a factor of 10, while the 'flash' 
is reduced by a factor of 50 because most of 
the higher-momentum pions are blocked by 
beam-line collimators. In the homogeneous B 
field of the storage ring, charged particles fol- 
low a circular path (slightly modified by the 
electric quadrupole field) which would cause 
them to strike the infiector after one revo- 
lution. Muon injection therefore requires an 
in-aperture magnetic pulse at 90" around the 
ring from the injection point (provided by the 
pulsed 'kicker') in order to center the muon 
orbits in the storage region. 

With either muon or pion injection, pos- 
itive (negative) muons are stored in the ring 
with spins initially polarized anti-parallel 
(parallel) to their momenta. In a magnetic 
field (no E-field) the spins precess relative 
to the muon momenta according to LUa = 
ujs — '^c = ^o,f^^. Here uja and iVc are the 
angular frequencies of spin rotation and mo- 
mentum rotation (or cyclotron angular fre- 
quency), respectively. Note that all of the 
muons precess at the same rate in a given 
field, regardless of their momenta. Two quan- 
tities must be measured with precision to de- 
termine a^: LUa and B, each being time av- 
eraged over the ensemble of muons. Actu- 
ally, instead of measuring B, we determine 
the frequency of precession of the free pro- 
ton, LUp, in the same average magnetic field 
as the muons via NMR measurements. The 
anomaly is given by Eq. 0, 



where R = A = ^ = 3.183 345 39(10) 

is the ratio of the muon and proton 
magnetic moments determined from other 
experimenttS'i. The analyses of < cJq > 
and < ujp > were independent, and further- 
more concealed offsets were maintained in 
each value so that no one could calculate 
prior to the completion of the analyses. 

The average trajectories of muons in the 
storage ring can only be known moderately 
well. Therefore the B field needs to be as 
uniform as possible to minimize the depen- 
dence of a given muon's precession rate on 
its exact trajectory in the storage ring. This 
prohibits the use of a gradient magnetic field 
to store (focus) the beam. E821 follows the 
CERN approach of using a quadrupole elec- 
trostatic field to provide the focusing. In the 
presence of the electric field, the precession is 
described by Eq. ^, 

= - — KB - (a,, - ^^)/3 X E] (3) 

The "magic" 7 « 29.3, or w 3.094 
GeV/c, is chosen so that — ~ 0, min- 
imizing the effect of E onuja- Because not all 
stored muons have the exact magic momen- 
tum, a small correction (« 0.6 ppm) must be 
applied to the final value for LUa- Equation 
3 is strictly valid only for muon motion per- 
pendicular to B; the up-down motion associ- 
ated with vertical betatron oscillations leads 
to another small correction of w 0.2 ppm, the 
so-called "pitch correction" . 

The practical limit to the strength of a 
ferric field with the required homogeneity is 
w 1.5T; E821 chose B = 1.45T, leading to 
a ring radius of 7.112m. The storage ring 
aperture radius is 4.5 cm, giving a « ±0.6% 
(« ±0.4%) base-to-base range in stored mo- 
menta for pion (muon) injection. The cy- 
clotron period is Tc = -J- = — w 149.2ns, the 



For Publisher's use 



llion e*{E >2GeV, 1999 dala) 



20 JO 40 bO oO 70 



Figure 12. Spectrum of number of positrons versus 
time, from the 1999 data sample. There are a total 
of 1 billion e+ above 2 GeV. 



precession period is Ta ~ 4.365 /xs and the 
dilated muon lifetime is r = 7x0 ~ 64.38 /iS. 
Decays are typically measured for at least ten 
muon lifetimes, or about 4000 cyclotron and 
150 precession periods. A log plot of the 1999 
data set, folded into 100/xs periods, is shown 
in Fig. |l|. 

The error on from the combined 
and data sets from the CERNEI g-2 ex- 
periment is 7 ppm with a 1.5 ppm system- 
atic error. By comparison, E821 must keep 
the systematic errors in 'Q{LOp) and uja to less 
than a few tenths of a ppm to approach its 
experimental precision goal of 0.35 ppm. 

3.1 The Magnet and the Determination 
ofujp 

The storage ring,@ Fig. ^ is a continuous 
C-magnet open to the inside. A cross-section 
view, Fig. [ij, shows its essential features. It 
contains more than 600 tons of magnet steel. 
Three superconducting coils, which provide 
exceptional B-field stability with time, are 
used to power the magnet. The entire magnet 
is wrapped in thermal insulation to reduce 




Figure 13. Overhead schematic view of the storage 
ring magnet. The detectors are distributed in uni- 
form intervals around the inside of the ring. The 
beam is brought in through a hole in the back of the 
magnet yoke, indicated by the solid line at 10 o'clock. 



gap changes due to temperature change. The 
storage region of 4.5 cm radius is defined by 
a series of circular collimators inside an evac- 
uated chamber. The pole gap is 18 cm high 
and 53 cm wide. 

Many shimming options were incorpo- 
rated in order to achieve the desired field uni- 
formity. The very high-quality steel of the 
pole tips is decoupled from the lower quality 
steel and the imperfections (including holes 
for cryogenic leads, etc.) of the yoke, by 
means of an air gap. Iron wedges in the gap 
can be moved radially to locally adjust the 
dipole field. The thickness and position of 
iron pole bumps can be adjusted to minimize 
quadrupole and sextupole fields. Thin sheets 
of iron were affixed to the pole tips to improve 
local uniformity. Current-carrying wires at- 
tached to circuit boards and mounted on the 
pole faces, with one set forming a closed loop 
covering 12*^ in azimuth and another set go- 
ing entirely around the ring, provide a final 
fine-tune of the dipole field. 

A continuous monitor of the B-field was 
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Figure 14. Cross section view of the storage ring mag- 
net. The means of shimming the magnet and mea- 
suring its field are indicated. 
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Figure 15. B-field contours across the storage region 
in E821, averaged in azimuth, for succeeding run cy- 
cles. The 1997 contours are 2 ppm, the rest 1 ppm. 



provided by 360 NMR probes placed in fixed 
positions around the ring, above and below 
the storage region. A subset of probes, those 
most highly correlated to the average B-field, 
provided feedback to the magnet power sup- 
ply to compensate for the slight field drift- 
ing which are mainly the result of ambient 
temperature changes. Two separate off-line 
analyses of the B-field used somewhat differ- 
ent combinations of probes to determine the 
average field as a function of time, with com- 
parable results. The B-field in the storage re- 
gion was mapped in 1 cm intervals every three 
to four days with 17 NMR probes mounted 
transversely on a movable cable-driven trol- 
ley. This was accomplished inside the vac- 
uum, with essentially no geometrical changes 
to the magnet or vacuum chamber configura- 
tion. The NMR probes on the trolley were 
calibrated against a standard spherical wa- 
ter NMR probe, which was normalized to the 
precession frequency of a free proton. The 
fixed probes, in the off-line analysis, tracked 
the trolley probes to better than 0.15 ppm 
over time. The steady improvement in the 
B-field provided by shimming is illustrated 
in Fig. The marked improvement from 
1999 to 2000 is attributable to the replace- 
ment of the inflector, which had a damaged 
superconducting fringe-field shield. 

The distribution of muons inside the stor- 
age region was determined from an analysis 
of the debunching of the beam as a function 
of time. At the time of injection, muons are 
localized in the ring with a full width at half 
maximum of about 120 degrees. As a re- 
sult, the time spectrum from a given detector 
at early times will contain oscillations with 
a period equal to the cyclotron period (the 
so-called "fast rotation" structure). Muons 
with high momenta have smaller cyclotron 
frequencies than those with low momenta, 
causing the bunches to spread out around 
the ring and the amplitude of the oscillations 
to diminish with time ( "debunching" lifetime 
~ 20 lis). The analysis of the debunching ver- 
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Table 3. Systematic errors in ujp for the 1999 data 
set. 



Source of errors 


Error (ppm) 


Inflector fringe field 


0.20 


Fixed probe calibration 


0.20 


Fixed probe interpolation 


0.15 


Trolley measurements 


0.10 


PL distribution 


0.12 


Absolute calibration 


0.05 


Otherst 


0.15 


Total Syst error on ujp 


0.4 



^ higher multipoles, trolley temperature sta- 
bility, kicker eddy currents. 



sus time gave the radius of curvature distri- 
bution of the muons, which in combination 
with simulations of the betatron motion of 
the muons produces the radial and vertical 
distribution of muons in the storage aperture. 
The distributions thus deduced were folded 
geometrically with the map of NMR frequen- 
cies to obtain < ujp >. Corroborating infor- 
mation on the horizontal and vertical distri- 
butions of muons, as well as information on 
the betatron motion, at early times, was pro- 
vided by scintillating fiber hodoscopes which 
could be inserted into the storage region. The 
hodoscopes were sufficiently thin that useful 
beam profile data could be taken for many 
tens of microseconds before the beam was de- 
graded. 

The final value for the average field is 
= 61 791 256 ± 25 Hz (0.4 ppm). The 
sources of systematic errors are given in Ta- 
ble ||. The improvements in the 2000 data set 
are the installation of a new inflector with far 
less fringe field, greatly reducing the first item 
in the Table, and better trolley calibrations. 



3.2 Determination of uj a 

The decay positrons from /i+ e'^VeU^ 
have energies in the range 3.1 GeV. 
In the muon rest frame, the higher energy 
positrons are preferentially emitted parallel 
to . When the muon spin is parallel to the 
muon momentum, there will be more high en- 
ergy muons in the lab frame than when the 
directions are anti-parallel. The number of 
positrons in the lab frame above a given en- 
ergy threshold Et versus time therefore oscil- 
lates at the precession frequency according to 
Eq. |, 

N{t) = Noe^i X 

(1+ AcOs{uJat + (l)a), (4) 

where each of A^o and A depend strongly on 
Et, while 0a depends slightly on Et- 

The positrons, generally having lower 
momenta than the muons, are swept by the 
B field to the inside of the storage ring, 
where they are intercepted by 24 scintillate 
ing fiber/lead electromagnetic calorimeter£3 
uniformly spaced around the ring. The typ- 
ical energy resolution of the calorimeters is 
~ E(GeV) • ^i'^*^^ ^ low energy positron 
arrives at the detectors more quickly after 
muon decay than a high energy positron (the 
average distance traveled is less), the actual 
measured times at the detectors relative to 
the muon decay time will depend slightly on 
energy, therefore 4'a depends slightly on Ef. 
4>a is highly correlated to Wa, implying the 
need to accurately calibrate the energy scale 
of the calorimeters. They were calibrated, on 
average, to better than 0.2%, using the ob- 
served energy spectra of decay positrons as 
a function of time. This corresponds to an 
0.02 ppm systematic error in Wo- The aver- 
age positron time measurement was stable to 
20 ps over any 200 /is time interval, as de- 
termined by a laser calibration system, also 
giving about an 0.1 ppm systematic error. 

Each calorimeter is equipped with 4 pho- 
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tomultiplier tubes, whose sum is sent to a 
waveform digitizer (WFD)l2I which samples 
the photomultipUer pulse height every 2.5 ns. 
Both the time of arrival and the energy of the 
positron are determined from the WFD infor- 
mation. 

For 1999, the analysis occurred in two 
steps. First, in the production step, WFD 
data were converted to positron energies and 
times. There were two separate productions 
of the data. They each developed indepen- 
dent algorithms to handle the WFD data, 
which eventually evolved to become similar. 
The second step involved performing ^ min- 
imization fits to the data in Fig. |l^ to obtain 
ijja- The parent distribution in Eq. ^ pro- 
vided a good fit to the 1998 data sample, 
with five variable parameters: iVo, r, A, uja 
and 0a. It did not however provide a good 
fit to the 1999 data, which has 15 times more 
positrons. It was necessary to account for 
small but noticeable effects from pulses over- 
lapping in time at high rates (pile-up), beta- 
tron motion of the stored muons, and muon 
losses. 

There were four independent analyses of 
the positron time spectra for the 1999 data 
set, two for each production. They used dif- 
ferent methods to handle these additional ef- 
fects. The time spectrum, with these addi- 
tional effects, can be described by a 14 pa- 
rameter function, Eq. ^ (not all parameters 
are necessarily variable in a fit): 



][t) = {iVoe--[l + Acoscjai-H(/)a] -Hp(i)} 
X h{t) X lit) (5) 

The pile-up term, p(t), with parameters 
rip, Ap, A0p, is given by Eq. I 



At times close to injection, the bunch- 
ing of the muons leads, in addition to the 
previously mentioned oscillations, to an en- 
hancement of the pileup. The oscillations are 
eliminated from the time spectra by using a 
bin width — Tc, adding to each arrival time 
a time uniformly randomized over ±-^, and 
summing all detectors around the ring. The 
pile-up enhancement is accounted for by the 
last term in Eq. ^; the constants are held fixed 
to values determined in separate pile-up stud- 
ies. 

Two of the analyses constructed a sim- 
ulated pileup time spectrum by combining 
single positron pulses, from data, into pile- 
up pulses. The pileup spectrum was then 
subtracted from the primary spectrum, thus 
eliminating the p(t) term from their fits. A 
third analysis varied rip and Ap, with A0p 
held fixed to the value determined by a fit 
to an artificial pileup spectrum. Fits made 
with A</)p variable produced a result consis- 
tent with other methods, but the statistical 
error on lOa was doubled as a consequence of 
the strong correlations of A(/)p to and 0q. 

The amplitude rip of pile-up was gener- 
ally less than 1% even at the earliest decay 
times, with an asymmetry small com- 
pared to the (g-2) asymmetry, A. The ar- 
tificial pileup spectrum gave the expected ^ 
lifetime, and when subtracted from the main 
spectrum, did a very good job of eliminating 
events above the maximum electron energy 
of 3.1 GeV, which apart from energy resolu- 
tion effects could only be due to pileup. The 
properties of the artificial pile-up spectrum 
matched very well with the results of multi- 
parameter pile-up fitting. 

The coherent betatron oscillation (CBO) 
term, b(t), with parameters Ab^ ^s, 0s, and 
TB^ is described by Eq. ItI, 



X (rip + Ay cos (uJoX + 0o + A0p)) 
X (l + apc"^'^^') (6) 



h{t) = 1 + As cos(wBt + 0b) 



(7) 



The need for b(t) is due to the effects 
of betatron motion of the muons combined 
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with the restricted aperture of the inflector. 
Muons are injected into the ring through the 
inflector, whose aperture was considerably 
smaller than the storage ring aperture be- 
cause of mechanical and geometrical limita- 
tions. This creates a muon beam with nar- 
row horizontal (1.8 cm) and vertical (5.6 cm) 
waists at the inflector exit at injection time. 
The more important horizontal waist case 
will be discussed here. In a perfectly uniform 
B field with no electric field, the muon trajec- 
tories, projected into the plane of the magnet, 
are circles. After the muons are kicked, the 
position of the horizontal waist would ideally 
be in the center of the storage region. The 
kick, however, was generally less than 100% 
of its optimum value. Therefore the average 
radius of muons at the narrow waist, at injec- 
tion time, was larger than the central storage 
region radius. At 180 degrees around the ring 
from the narrow waist, the muons are spread 
out to fill the ring aperture, and have an av- 
erage radius more nearly equal to the central 
ring radius. The acceptance of the electron 
calorimeters depends to a slight extent on the 
horizontal width and especially on the aver- 
age radial position of the muon beam. Thus 
the detector positron acceptance will be dif- 
ferent at the narrow waist compared to the 
opposite side of the ring. When we add elec- 
trostatic focusing, the position of the narrow 
waist (or focus) will move around the ring 
at the so-called CBO frequency, which is the 
cyclotron frequency minus the horizontal be- 
tatron frequency, fcBO = /c(l - \/l - n). 
For the field index n = -^j^^ = 0.137, 

/3Bo Or ' 

fcBO ~ 475 kHz. As a result, we get a small 
oscillation at « 475 kHz superimposed on the 
time spectrum. The amplitude of the CBO 
is typically a few tenths of a percent of the 
total number of counts, and Tf, is long, about 
100 fis. 

One analysis allowed all four CBO pa- 
rameters to vary. Two kept lub fixed to the 
value obtained from a Fourier transform of 
the residuals from a five-parameter fit. The 



frequencies from the Fourier transform and 
the fit were in good agreement. 

The muon loss term, Eq. |[ has two pa- 
rameters, and t^l. 



l{t) = 1-Fa^L -e ^^'^^ 



(8) 



Muon losses are thought to be caused 
by the slight drift of the orbits of muons 
whose trajectories bring them close to col- 
limators. The drift could be caused by the 
small non-uniformities in the E and B- fields, 
although the exact mechanism is not known. 
Indeed, when the beam is 'scraped' for about 
15 IIS right after injection, by temporarily 
displacing the stored muon beam several mil- 
limeters vertically and horizontally in order 
to force the loss of muons with trajectories 
close to the collimators. The rate of muon 
loss is markedly reduced after the scraping 
is turned off, compared to the no scraping 
case. After scraping, losses are generally less 
than ^1% at early decay times, with the rate 
of losses decreasing with a short lifetime of 
Tfj_L ~ 20/is. In addition there was a roughly 
constant loss rate of « 0.1% per muon life- 
time at all decay times, as determined by 
comparing the measured decay rate with that 
expected from special relativity. The two 
analyses which allowed both parameters to 
vary obtained the same loss lifetime as was 
observed in a third analysis which held r^^ 
fixed to the value determined using separate 
muon loss detectors. 

It is important to realize that only two of 
the 14 parameters in Eq. ^ have a large corre- 
lation to uJa'- 4>a and A(/)p. The latter param- 
eter received a great deal of study during the 
pile-up analysis. The former introduces no 
systematic error, but increases the statistical 
error by a factor 

The fourth analysis was the novel "ratio 
fit" . After pile-up subtraction, events were 
randomly assigned to four separate time his- 
tograms. Then, the ratio in Eq. S was formed: 
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Nl 



N2 



Table 4. Results of the four analyses for oJa- -R is 
defined by oja = 27r/o(l — -Rx 10~®). /o is the nominal 
precession frequency. 



# Par. 


xVdof 


R(ppm) 


13 


1.012 ±0.023 


143.24 ± 1.24 


10 


1.005 ±0.023 


143.08 ± 1.24 


9 


1.016 ±0.015 


143.30 ± 1.23 


3 


0.986 ±0.025 


143.37 ± 1.28 


Avg. 




143.17 ± 1.24 



Table 5. Systematic errors in uia for the 1999 
set. 



data 



Source of errors 


Error (ppm) 


Pile- up 


0.13 


AGS background 


0.10 


Lost niuons 


0.10 


Timing shifts 


0.10 


E field, pitch 


0.08 


Binning, fit procedure 


0.07 


Debunching 


0.04 


Gain changes 


0.02 


Total Syst error on uja 


0.3 



N+ + + Nl + N2 

= AcosiiUat + M + i-^f (9) 
16r^ 

where N+ = Ni{t + N- = N2{t - 
AfO = iV3(t), and = Ni{t). The muon 
lifetime cancels, and r{t) is sufficiently insen- 
sitive to the CBO and the muon losses that 
these effects can be neglected in the fit. The 
insensitivity to the GBO is a consequence of 
— being not so different from We arrive 
at a three parameter fit which has different 
responses to systematic errors compared to 
the conventional multi-parameter fits. 

The results of the four analyses are given 
in Table 0. All of the results are well within 
the bounds expected for correlated data sets. 
The final value for tOa is the average of 
these results, ^ = 229072.8 ± 0.3 Hz(1.3 
ppm), after a correction of +0.81 ± 0.08ppm 
for the effects of the electric field in Eq. |^ 
and for vertical betatron oscillations ( "pitch" 
correction) .□ Note that the data are not sen- 
sitive to the sign of uJa] this however is well- 
determined from many other experimental 
measurements, and is implicit in the value 
of A which we use to extract the value of 
a^. The sources of systematic error in oja are 
given m Table |. The AGS background is the 
result of unwanted particle injection into the 
ring after the initial injection. The improve- 
ments for the 2000 data run were the addition 
of a sweeper magnet in the beam-line to elim- 
inate errors due to the AGS background, and 
an increase in the number of lost muon detec- 
tors in order to reduce the muon loss error. 



4 Conclusions and Outlook 

The muon anomalous magnetic moment can 
be both measured and calculated (within the 
Standard Model) to a high precision, and 
given its high sensitivity to new physics, its 
measurement affords an exceptional oppor- 
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tunity to probe for new physics beyond the 
Standard Model. 

The new world average value of a+ shows 
a 2.2 ± 1.4 ppm difference from the Czarnecki 
and Marcianoi theory compilation, after the 
sign of the LOL pole term is corrected. Many 
theoretical ideas have been put forward to ex- 
plain any difference, including supersymme- 
try, leptoquarks, muon substructure, etc. It 
could of course also be explained by a statisti- 
cal fluctuation, an error in the experiment, or 
an error in the Standard Model calculation. 

All aspects of the theoretical calculation 
of are being heavily scrutinized. New high 
quality e+e~ data from VEPP-2M and Bei- 
jing as well as r decay data from LEP, are be- 
ing analyzed now, and should have an impact 
on a*/* in the next few months. Longer term, 
Novosibirsk and Beijing have upgrade plans, 
and DA<i>NE (and perhaps the B-factories 
and Cornell) have plans to use the radia- 
tive return process to measure e+e~ cross 
sections. Further calculations of the light- 
by-light term are being considered by several 
groups. One can reasonably expect a con- 
tinued steady improvement in the error and 
reliability of a^''. 

In E821, analysis is under way on the ~ 
4 billion positrons (/i^) from the 2000 run 
(about four times larger than the 1999 data 
set) and on the ~ 3 billion electrons (/i~) 
from the 2001 run. Systematic errors are 
expected to be reduced for both data sets. 
Once these data are analyzed, it should be 
possible to make a more definitive statement 
concerning whether the measured anomaly 
agrees with theory. Comparison of the /i"*" 
and fi~ anomalies is a test of both the sys- 
tematic errors in E821 and also CPT invari- 
ance. With another data run, E821 expects 
to achieve 0.3 ppm statistical error on a'ifP 
and an estimated 0.3 ppm systematic error, 
not far from the original goal of 0.35 ppm 
overall error. 

It is interesting to note that any new 
physics affecting may also lead to a non- 



zero permanent electric dipole moment for 
the muon, through its CP violating part. 
Assuming that the CP violating phase for 
new physics (I)cp ~ Ij then dimensional ar- 
guments, along with the observed value for 



give 



10 



'22 



Even if exper- 



iment and theory were to agree, the muon 
EDM is interesting in its own right: it is 
the only currently accessible EDM from a 
second generation particle. Comparing with 
the electron, the current limit on the electron 
dipole moment is « 4 x 10~^^e — cm. If the 
EDM scales by the first power of the mass, 
then a lO^^^e — cm muon measurement is 
competitive with that of the electron. There 
are speculations, however, that the electron 
EDM could be small due to an accidental can- 
cellation which may not apply to the muon. 
Or, if the scaling is with the square of the 
mass or higher, the muon then becomes more 
sensitive than the electron to new physics.Ej 
There are a number of models which predict 
in the range 10~^^e — cm to 10~^''e — cm.LJ 
The presence of an EDM adds the term 



er] 



{(3xB + E) 



(10) 



2mc 

to Eq. H, where the EDM is given by = 
^(imc)- ^^"^ (g"2) experiment, the ef- 
fect of the dominant f3 x B term is to tip 
the precession vector radially by an angle 
B = tan^^ TT^. This causes an increase in 



the precession frequency to « uja\/^ + rf. 
It also causes an oscillation about zero of the 
average vertical component of the positron 
momenta, which can be observed as an oscil- 
lation, with frequency lu, in the average ver- 
tical position of positrons on the face of the 
calorimeters. 

In the unlikely event that all of Aa^ 
can be attributed to a muon EDM, then 
d^. = (2.3 ± 0.7) X 10-13. The CERN (g- 
2)eI experiment has set the best limit on the 
muon EDM so far, < 1 x I0~^^e — cm, de- 
duced from limits on the vertical oscillations. 
While this value for is larger than any the- 
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ory predicts, it is nevertheless not ruled out 
by CERN limit. E821 expects to reduce the 
limit by about a factor of five from an im- 
proved measurement of the vertical oscilla- 
tions. 

A dedicated experiment to measure the 
muon EDM to the 10~^"'e — cm level is cur- 
rently being develDped at Brookhaven Na- 
tional Laboratory .cll It would use a new tech- 
nique where a muon momentum and an ap- 
plied electric field would be selected so that 
the second term cancels the first term in 
Eq. 1^. One is only left with the motion in 
a vertical plane described by Eq. In this 
technique, there is a very large enhancement 
of the EDM signal relative to the "noise" over 
the technique used in the g-2 experiments. It 
is planned to mount this experiment over the 
next few years. 
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